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ABSTRACT 


A  computer  routine  composed  of  two  programs  was  developed  to 
process  and  analyze  aircraft  reconnaissance  reports  for  the  purpose  of 
diagnosing  the  structure  of  the  tropical  cyclones  of  the  eastern  North 
Pacific  Ocean  (EASTROPAC) .     The  reports  processed  by  the  routine  are 
the  Reconnaissance,  Detailed  Vortex/Center,  and  the  Dropsonde  Reports. 

Data  from  the  year  1971  were  used  as  a  prototype  set  for  the  develop- 
ment of  the  computer  routine.     The  reconnaissance  observations  were  de- 
coded by  the  first  program  and  analyses  were  drawn  by  the  second  program 
at  various  levels  between  1000-mb  and  200-mb  during  the  storm  and 
hurricane  stages  of  the  tropical  cyclones.     The  data  proved  to  be  adequate 
for  the  purpose  of  proving  the  workability  of  the  computer  routine,  but 
insufficient  for  anything  but  a  preliminary  analysis  of  EASTROPAC  cyclone 
structure. 
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I .      BACKGROUND  AND  OBJECTIVES 

Since  the  launching  of  the  first  operational  weather  satellite  in  1964, 
interest  in  the  cyclones  of  the  eastern  tropical  North  Pacific  Ocean 
(EASTROPAC)  has  been  on  the  increase.     Prior  to  that  time,  a  tropical 
vortex  was  usually  well  developed  before  its  existence  was  even  suspected, 
and  almost  all  information  pertaining  to  these  cyclones  was  obtained  from 
ships  which  had  inadvertently  wandered  into  their  paths.     The  data  were, 
at  best,   sparse  and  unprofessionally  observed,   being  useful  mainly  for 
the  positioning  of  the  cyclone. 

With  the  advent  of  timely  satellite  information,  cyclones  could  be 
detected  in  the  early  stages  of  development,  and  military  reconnaissance 
flights  could  be  initiated  to  provide  not  only  confirmation  of  position,  but 
also  real-time  measurement  of  the  meteorological  parameters  in  the 
vicinity  of  the  cyclone.     To  the  knowledge  of  this  researcher,  use  of 
these  data  have  not  been  made  for  purposes  of  EASTROPAC  research,  and, 
in  fact,  no  comprehensive  study  of  the  structure  of  EASTROPAC  cyclones 
has  been  undertaken  to  date. 

Simpson  and  Frank  (1970)  suggested  that  the  EASTROPAC  area  be 
treated  merely  as  an  extension  of  the  North  Atlantic/Carribbean  area, 
which  implies  a  similarity  of  the  structure  of  EASTROPAC  cyclones  with 
those  of  the  Atlantic;  however,  Hansen  (1972)  has  shown  through  clima- 
tology that  there  are  significant  differences  in  size,  frequency,  and 
behavior  between  EASTROPAC  cyclones  and  those  of  the  Atlantic,  as  well 
as  the  western  Pacific  Ocean  areas. 


Serra  (1971)  states  that  inadequate  data  exist  at  present  for  the  study 
of  these  cyclones  and  suggests  that,  as  an  interim  measure,  the  knowl- 
edge of  the  hurricanes  of  the  Atlantic  area  be  applied  to  studies  of  the 
EASTROPAC  cyclones. 

This  researcher  feels  that  the  information  collected  by  military  recon- 
naissance flights  since  1964  (perhaps  in  combination  with  ship  and 
satellite  data)  should  provide  an  initial,   if  not  adequate,  data  base  for 
research  into  the  structure  of  the  EASTROPAC  cyclones  without  undue  re- 
liance on  cyclone  models  developed  for  other  ocean  areas.     The  processing 
and  presentation  of  the  appropriate  aircraft  reconnaissance  data  for  the 
aforementioned  purpose  is  the  major  objective  of  this  study. 


II.      THE  NATURE  OF  THE  PROBLEM 

For  the  purpose  of  this  study,  the  term  "structure"  ,  as  applied  to 
EASTROPAC  cyclones,  will  be  limited  to  the  horizontal  and  vertical  dis- 
tribution of  wind,  dry-bulb  and  dew-point  temperatures,  and  the  deviation 
of  the  isobaric  heights  from  the  United  States  Standard  Atmosphere  values. 
These  particular  parameters  were  chosen  because  of  their  availability  and 
their  relationship  to  the  important  dynamic  and  thermodynamic  processes 
of  the  tropical  cyclone.     Some  or  all  of  these  parameters  are  included  in 
each  of  the  three  major  types  of  meteorological  reports  submitted  by 
military  aircraft  on  cyclone  reconnaissance  missions,  namely  the  Recon- 
naissance Report  (RECCO)  ,  the  Detailed  Vortex/Center  Report  (DETAILED 
EYE) ,  and  the  Dropsonde  Report  (DROPSONDE) . 

The  data  from  the  reports  just  enumerated  are  available  in  a  coded 
format  from  the  meteorological  worksheets  carried  by  the  aircraft  on  re- 
connaissance missions.     These  data  exist  for  varying  altitudes  (flight 
levels)  and  for  all  stages  of  development.     Before  these  raw  data  can  be 
used,  they  must  be  decoded  and  reduced  to  a  homogeneous  form.     Then 
they  must  be  properly  assembled  and  analyses  drawn  before  any  determina- 
tion can  be  made  of  the  structure  of  the  EASTROPAC  cyclones.     In  view  of 
the  large  variation  of  the  data  in  both  time  and  space,  the  manner  in  which 
they  are  to  be  assembled  in  order  to  give  the  most  meaningful  result  be- 
comes a  major  consideration.     The  most  straightforward  approach  might  be 
to  assemble  all  available  data  taken  at  (or  near)  a  given  level  into  a 


single  composite  cyclone,  or,  depending  on  data  density,  the  same  data 
might  be  assembled  into  a  composite  cyclone  at  each  of  the  stages  of 
development.     The  list  could  go  on  and  on,  and  each  possibility  has  its 
own  particular  merits. 

Because  of  the  sheer  volume  of  data  available  and  the  multiplicity  of 
ways  in  which  it  could  be  assembled,  the  approach  which  best  facilitates 
handling  the  data,  and  hence  relates  to  obtaining  the  most  meaningful  re- 
sults, is  through  the  use  of  the  digital  computer.     (The  advantages  of 
using  a  modern  high  speed  computer  for  handling  large  quantities  of  data 
in  a  diverse  manner  are  well  known  and  will  not  be  discussed  here.) 
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III.      THE  COMPUTER  ROUTINE 

While  a  single  program  might  be  written  which  would  satisfy  the  objec- 
tives of  this  study,  that  program  would  require  vast  amounts  of  core  stor- 
age and  central  processing  unit  (CPU)  time  for  each  run.    Considerable 
savings  in  CPU  time  could  be  gained  if  the  job  were  divided  so  that  the 
data  were  decoded  only  once  rather  than  each  time  a  different  investiga- 
tion is  made.    For  these  reasons  it  was  determined  that,  in  terms  of  over- 
all efficiency,  a  computer  routine  made  up  of  two  separate  programs  (here- 
inafter called  PROGRAM  ONE  and  PROGRAM  TWO)  would  best  satisfy  the 
objectives  of  this  study. 

A.      PROGRAM  ONE 

PROGRAM  ONE  consists  of  a  main  program  and  six  subroutines.     It 
performs  three  major  functions:     i)  the  decoding  of  the  RECCO,   DETAILED 
EYE,  and  DROPSONDE  reports,  ii)  the  transformation  of  geographic  coor- 
dinates into  moving  natural  coordinates,  and  iii)  the  punching  of  data 
cards  for  PROGRAM  TWO.     It  takes  as  input  punched  cards  produced  from 
the  aircraft  reconnaissance  worksheets  and  produces  as  output  both 
printed  tables  and  punched  cards. 

The  program  first  reads  the  six-hourly  cyclone  center  (BEST  TRACK) 
positions  obtained  from  the  applicable  Fleet  Weather  Central  Guam, 
Annual  Typhoon  Report  and  the  associated  maximum  sustained  wind  speeds 
taken  from  the  official  forecast  bulletins  issued  by  the  cognizant  Fleet 
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Weather  Central.*    This  information  is  then  used  to  compute  the  course 
and  speed  of  the  cyclone  for  each  six-hour  segment  of  the  track,  along 
with  the  rate  of  increase  or  decrease  of  the  maximum  wind  speed  for  each 
six-hour  interval.     This  information  is  printed  in  the  form  of  a  table, 
(see  Table  I)  and  is  stored  for  use  in  the  transformation  of  coordinates 
(see  below)  . 

The  three  types  of  coded  reports  are  then  read,  and  the  RECCO  Report 
decoded  in  accordance  with  procedures  found  on  the  RECCO  worksheet, 
the  DETAILED  EYE  Report  decoded  in  accordance  with  procedures  found  in 
the  "National  Hurricane  Operations  Plan"  [6],  and  the  DROPSONDE  Report 
decoded  in  accordance  with  procedures  found  in  the  United  States  Air 
Force  Ninth  Weather  Reconnaissance  Wing  Manual  105-1,  Vol.  I  [5]. 
Realizing  that  future  research  might  utilize  other  data  sources,   (satellite, 
ship  reports,  etc.)  the  program  was  written  with  enough  flexibility  so 
that  up  to  95  additional  types  of  reports  could  be  processed  by  making 
only  minor  changes  to  the  program.     During  the  decoding  process  gross- 
error  and  missing-data  checks  are  made  and  such  data  identified  to  the 
program. 

The  effect  of  the  geographical  position  of  the  cyclone  system  is  then 
removed  through  a  transformation  of  coordinates.     The  geographic  coor- 
dinates (latitude  and  longitude)  of  each  report  are  transformed  into  a 


*    Fleet  Weather  Central  (now  Fleet  Weather  Facility)  ,  Alameda, 
California  for  EASTROPAC  cyclones  east  of  140W;  Fleet  Weather  Central, 
Pearl  Harbor,  Hawaii,  for  EASTROPAC  cyclones  west  of  140W. 
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natural  coordinate  system  whose  origin  is  at  the  computed  position  of  the 
center  of  the  vortex  at  the  time  of  the  report.     This  position  is  computed 
from  BEST-TRACK  data,  and  the  transformation  is  accomplished  in  the 
following  Manner:    The  BEST  TRACK  is  searched  to  locate  the  proper  six- 
hour  segment,  and  the  center  of  the  vortex  is  located  along  that  segment 
at  a  distance  equal  to  the  difference  between  the  time  of  the  last  BEST- 
TRACK  position  and  the  time  of  the  report  multipled  by  the  average  speed 
of  movement  of  the  center  along  that  segment.    With  the  location  of  the 
origin  thus  fixed,  the  position  of  the  report  is  converted  into  polar  coor- 
dinates (r,  9)  ,  where  r  is  the  radial  distance  in  nautical  miles  from  the 
center  of  the  vortex  to  the  position  of  the  report  and  8  is  the  angle  meas- 
ured clockwise  from  true  north  to  the  radial  line  joining  the  center  of  the 
vortex  and  the  position  of  the  report.     See  Figure  1.     The  entire  polar 
coordinate  system  is  then  rotated  through  the  angle  0,  the  smaller  angle 
between  the  direction  of  vortex  movement  and  true  north,   so  that  the 
direction  toward  which  the  cyclone  is  moving  coincides  with  true  north. 
See  Figure  2.  At  this  point  the  position  of  the  report  relative  to  the 

center  is  then  reconverted  into  a  Cartesian  Coordinate  System  in  which 
the  distance  along  both  axes  is  measured  in  nautical  miles.    (This  interim 
conversion  to  polar  coordinates  was  made  only  to  facilitate  the  rotation 
of  the  axes,   since  rotation  is  accomplished  more  easily  in  a  polar  coor- 
dinate system  than  in  a  Cartesian  Coordinate  System.)     In  terms  of  this 
new  Cartesian  Coordinate  System,  the  direction  toward  which  the  cyclone 
is  moving  represents  positive  values  of  Y  and  directions  to  the  right  of  the 
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Figure  1.     Polar  and  geographic  coordinate  systems 


Figure  2.    Cyclone-centered  natural  coordinate 
system 
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track  represent  positive  values  of  X.     In  current  cyclone  nomenclature, 
the  right-front,  left -front,  left -rear,  and  right-rear  quadrants  of  the  cyclone 
correspond,  respectively,  to  the  first  (I),   second  (II),  third  (III),  and 
fourth  (IV)    quadrants  of  the  new  Cartesian  Coordinate  System.     Therefore, 
in  Figure  3,  the  schematically  represented  report  lays  in  the  right-front/I 
quadrant.    The  net  effect  of  this  transformation  is  the  creation  of  a  moving 
natural  coordinate  system  whose  origin  is  at  the  center  of  the  vortex,  and 
whose  speed  and  direction  of  movement  is  that  of  the  cyclone  itself. 

The  use  of  BEST-TRACK  data  in  the  transformation  procedure  is  be- 
lieved to  be  quite  accurate,   since  in  every  case  where  a  comparison  could 
be  made  between  a  vortex  center  position  computed  from  BEST-TRACK  data 
and  an  actual  aircraft  fix  of  the  center,  the  difference  between  the  com- 
puted position  and  the  reported  position  was  less  than  the  aircraft's 
estimated  error  in  navigational  position.     (In  the  DETAILED  EYE  Report  the 
location  of  the  center  is  reported  along  with  an  estimate  of  the  accuracy  of 
the  navigational  fix,  usually  ranging  from  2-50  miles.)    For  this  reason, 
and  to  preserve  continuity  in  the  calculations,  computed  positions  are 
used  in  the  transformation  of  coordinates  even  though  an  actual  fix  might 
be  available  for  the  time  of  a  report. 

Once  the  decoding  and  transformation  of  coordinates  is  completed, 
the  program  then  produces  an  output  of  one  or  more  punched  cards,  ac- 
cording to  the  type  of  report.     Each  RECCO  Report  produces  one  card,  each 
DETAILED  EYE  Report  produces  up  to  four  cards,  and  each  DROPSONDE 
Report  produces  up  to  seven  cards.     These  cards  contain  the  decoded  data, 
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with  all  positions  given  in  terms  of  the  natural  coordinate  system,  and 
are  used  as  data  for  PROGRAM  TWO.    Additionally,  a  table  is  printed 
containing  the  decoded  information  (see  Table  II)  .     In  this  table  the  posi- 
tions are  given  in  the  geographic  coordinate  system  as  well  as  in  the 
natural  coordinate  system. 

B.      PROGRAM  TWO 

PROGRAM  TWO  consists  of  a  main  program  and  three  subroutines.     It 
also  performs  three  major  functions:    i)  storing  the  data,   ii)  assembling 
the  data,  and  iii)  analyzing  the  data.     It  takes  as  input  the  decoded  data 
cards  produced  by  PROGRAM  ONE  and  produces  as  output  both  printed 
tables  and  analyses  drawn  on  the  CALCOMP  plotter. 

The  program  reads  the  data  and  stores  them  in  a  series  of  vectors, 
with  a  separate  vector  for  each  type  of  information  (e.g.  ,  one  for  wind 
speed,  one  for  wind  direction,  etc.).    At  the  option  of  the  user,  a  table 
may  be  produced  at  this  point  containing  an  "echo  check"  of  the  data 
cards  as  they  are  read. 

Next,  a  check  is  made  to  determine  whether  or  not  flight  level  analyses 
are  desired  and,  if  so,  for  which  flight  level.     If  no  flight  level  plots  are 
requested,  the  program  checks  to  see  whether  D-value  analyses  are  de- 
sired as  described  below.     If  flight  level  plots  have  been  requested,  a 
search  is  then  begun  through  the  vector  containing  the  flight-level  data. 
When  a  flight  level  is  found  which  lays  within  plus  or  minus  500  meters 
of  the  requested  level  another  check  is  made  to  determine  whether  or  not 
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the  point  lays  within  the  specified  distance  of  the  center  (see  below)  .     If 
the  point  is  found  to  lay  near  enough  to  the  center,    (i.e.  within  the  grid) 
further  tests  may  be  made,  at  the  option  of  the  user,   for  stage  of  develop- 
ment, year,  month,  etc.     This  optional  test  allows  the  data  to  be  assem- 
bled in  an  almost  limitless  number  of  combinations,  according  to  the  type 
of  research  being  conducted.     Once  all  tests  have  been  successfully  com- 
pleted, the  wind  speed,  wind  direction,  dry-bulb  and  dew-point  tempera- 
tures, and  the  positional  information  are  extracted  and  set  into  new  storage 
arrays. 

The  same  tests  are  conducted  for  each  point  in  the  flight -level  vector 
until  all  points  which  meet  the  specified  criteria  have  been  identified  and 
shifted  into  the  new  storage  array.    At  this  point,  the  assembly  of  the 
data  is  complete.     The  data  are  stored  separately  and  are  ready  to  be 
passed  to  the  subroutines  for  analysis. 

The  actual  analysis  is  performed  in  two  steps.     In  the  first  step  the 
data  points,  whose  positional  information  is  given  in  terms  of  Cartesian 
Coordinates  (X,Y),  are  passed  to  subroutine  GRID,  which  is  an  integral 
part  of  PROGRAM  TWO.     This  subroutine  takes  these  randomly  placed  data 
points  and  interpolates  them  into  a  21  x  21  grid  which  is  centered  over  the 
center  of  the  vortex.     This  grid  is  then  passed  to  the  library  subroutine 
CONTUR  which  draws  the  actual  analysis  using  the  CALCOMP  plotter. 

A  check  is  then  made  to  determine  whether  or  not  D-value  analyses 
are  desired,  and  a  similar  process  is  followed  in  selecting  these  points 
and  performing  the  analysis.    In  this  case,  the  selection  criterion  specifies 
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that  the  report  point  lay  exactly  at  the  required  pressure  level  rather  than 
within  some  specified  pressure  range  of  that  level. 

During  a  single  run  of  PROGRAM  TWO,  any  combination  of  the  follow- 
ing analyses  may  be  generated  at  the  option  of  the  user: 

a.  At  a  specified  flight  level: 

(1)  isotach  analysis  (10  kt  interval) 

(2)  isotherm  analysis  (5  C  interval) 

(3)  isodrosotherm  analysis  (5  C  interval) 

b.  At  any  two  isobaric  surfaces: 

(1)     D-value  analysis  (50  m  interval) 
Again,  at  the  option  of  the  user,  a  table  may  be  generated  showing  the 
data  points  selected  for  use  in  each  of  the  analyses. 

The  isotach,   isotherm,  and  isodrosotherm  analyses  are  performed  in 
the  vicinity  of  a  particular  altitude  rather  than  at  a  specific  pressure  level, 
since  current  reporting  procedures  specify  that  wind  speed  and  dry-bulb 
and  dew-point  temperatures  be  reported  by  the  aircraft  at  flight  level  (in 
meters)  rather  than  as  values  extrapolated  to  the  nearest  standard  pressure 
surface . 
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IV.      THE  PROCEDURE 

The  year  1971  was  arbitrarily  chosen  as  a  prototype  year  for  the  study. 
The  Naval  Weather  Service  Environmental  Detachment,  National  Climatic 
Center,  Asheville,  North  Carolina,  provided  copies  of  the  worksheets  sub- 
mitted by  United  States  Air  Force  aircraft  on  cyclone  reconnaissance  flights 
during  that  year.     These  sheets  yielded  471  RECCO  Reports,   28  DETAILED 
EYE  Reports,  and  59  DROPSONE  Reports,  for  a  total  of  558  reconnaissance 
reports  of  all  types.    All  available  information  leads  this  researcher  to 
believe  that  this  represents  approximately  95%  of  all  cyclone  data  recorded 
by  EASTROPAC  reconnaissance  aircraft  in  1971.     The  BEST-TRACK  informa- 
tion was  obtained  from  the  Fleet  Weather  Central  Guam  Annual  Typhoon 
Report  1971  [7],  and  the  maximum  sustained  wind  information  was  taken 
from  the  official  forecast  bulletins  issued  by  Fleet  Weather  Central  (now 
Fleet  Weather  Facility)  Alameda,  California. 

PROGRAM  ONE  was  then  tested  using  a  small  portion  of  actual  data 
from  the  year  1971 .     In  order  to  minimize  the  time  required  to  test  the 
analysis  scheme,  a  contrived  data  set  was  used  to  test  PROGRAM  TWO. 
The  results  of  this  test  are  shown  in  Figures  3  and  4.     Figure  3  is  a  sam- 
ple isotach  analysis  produced  from  174  data  points.    A  hand-drawn  (dashed 
line)  analysis  has  been  superimposed  on  the  computer-drawn  (solid  line) 
analysis  for  comparison.     Figure  4  is  a  sample  D-value  analysis  produced 
by  PROGRAM  TWO  from  59  data  points.     For  comparison,  the  D-values  at 
each  of  the  data  points  have  been  superimposed  on  the  computer  drawn 
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Figure  3.     Sample  isotach  analysis  drawn  from  174  contrived  data  points 
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Figure  4.     Sample  D-value  analysis  drawn  from  59  contrived  data  points 
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analysis.     These  figures  serve  to  illustrate  the  accuracy  with  which  the 
analyses  can  be  performed  given  a  sufficiently  dense  data  distribution. 

Once  the  programs  were  tested  and  the  parameters  to  be  investigated 
had  been  chosen,  the  information  contained  on  the  worksheets  had  to  be 
reduced  to  a  form  suitable  for  computer  processing.     A  series  of  formats 
was  devised,  and  each  type  of  coded  report  (in  raw  form)  was  reduced  to 
one  or  more  punched  cards  for  input  to  PROGRAM  ONE.     The  data  were 
manually  verified  from  the  printed  tables,  and  a  processed  data  set  was 
punched  by  PROGRAM  ONE. 

With  this  data  set  as  input  for  PROGRAM  TWO,   investigations  were 
made  into  the  structure  of  the  1971  cyclones  during  both  tropical-storm 
and  hurricane  stages.     (The  tropical-depression  stage  was  not  investiga- 
ted and  is  left  as  a  topic  for  future  research.)    In  each  investigation  a 
mesh  length  of  20  nautical  miles  was  used  in  both  the  X  and  Y  directions. 
Thus,  the  area  over  which  the  analyses  were  performed  was  a  square, 
400  nautical  miles  on  a  side,  centered  over  the  computed  (BEST-TRACK) 
position  of  the  composite  cyclone.     Selection  of  this  area  was  based  on 
the  assumption  that,  for  the  parameters  used,  the  contribution  to  the  de- 
lineation of  the  structure  of  the  cyclone  by  reports  which  lay  at  a  distance 
of  greater  than  200  nautical  miles  from  the  center  of  the  vortex  would  be 
negligible.     Isotach,   isotherm,  and  isodrosotherm  analyses  were  made  at 
flight  levels  of  3000  meters  (near  700  mb)  ,   5600  meters  (near  500  mb)  , 
7200  meters  (near  400  mb)  ,  and  9200  meters  (near  300  mb) .     D-value 
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analyses  were  made  at  the  1000-mb,   850-mb,   700-mb,   500-mb,  400-mb, 
300-mb,  and  200-mb  pressure  levels,  using  the  United  States  Standard 
Atmosphere  heights  of  each  of  the  surfaces  as  the  base  value  (i.e.  D=0) . 

Before  the  gridding  of  the  data  points,  the  fields  were  initialized  to 
the  following  values: 


Isotachs 

15  kt 

Isotherms 

0  C 

Isodrosotherms 

0  C 

D-values 

0  m 

This  initialization  scheme  applies  very  little  bias  to  the  data,  allowing 
the  computer  to  generate  analyses  which  are  largely  determined  by  the 
actual  data  while  conserving  CPU  time  in  the  gridding. 

The  results  of  these  investigations  appear  in  the  following  section. 
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V.      THE  RESULTS 

After  processing  by  PROGRAM  ONE,  the  558  reconnaissance  reports 
from  the  year  1971  yielded  898  data  points  to  be  used  in  the  analyses. 
The  following  table  contains  a  summary  of  the  number  of  points  which 
were  found  to  be  valid  for  each  analysis  within  the  previously  stated 
selection  criteria  for  each  level  investigated  in  the  study. 


Table  HI.     Summary  of  1971  Data 


Stage 


3000  meters 
Isotach 
Isotherm 
Isodrosotherm 

5600  meters 
Isotach 
Isotherm 
Isodrosotherm 

7200  meters 
Isotach 
Isotherm 
Isodrosotherm 

9200  meters 
Isotach 
Isotherm 
Isodrosotherm 

D-value  Analyses 

1000  mb 
850  mb 
700  mb 
500  mb 
400  mb 
300  mb 
200  mb 


Tropical  Storm 

23 
43 
27 


4* 
11* 

6 


20 
21 
46 
14 

9 
20 

0 


Hurricane 

21 
38 
23* 


8 

14 

6 


4 
6 
3 


1 
18* 
46 
16 

7 
11 

0 


The  analyses  marked  with  an  asterisk  ('*)  appear  as  Figures  5-8 
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Several  prototype  analyses  are  presented  (Figures  5-8)  ,  more  for  evi- 
dence of  a  viable  analysis  routine  than  for  diagnosis  of  tropical-cyclone 
structure.     The  following  is  a  discussion  of  certain  features  of  these 
analyses. 

The  isodrosotherm  analysis  in  Figure  5  illustrates  one  type  of  problem 
which  may  be  encountered  in  the  use  of  a  strictly  objective  analysis 
scheme  over  a  sparse  data  area.     The  very  dry  cell  analyzed  in  quadrant  I, 
near  the  center  of  the  figure,   is  believed  to  have  been  generated  almost 
entirely  by  the  -34  C  report  laying  in  the  second  quadrant  approximately 
20  nautical  miles  from  the  center.     The  value  at  this  single  point  is  more 
than  40  C  lower  than  the  value  at  any  other  point  in  the  first  or  second 
quadrant  except  those  which  lay  in  close  proximity  to  it.     This  fact, 
coupled  with  its  position  relative  to  the  other  points  in  these  quadrants, 
caused  the  analysis  scheme  to  generate  the  closed  dry  cell  in  order  to 
obtain  the  proper  gradient.     The  -32  C  report  (second  quadrant,  at  approxi- 
mately 15  nautical  miles  from  the  center)  did  not  cause  the  same  type  of 
cell  to  be  generated  because  of  the  absence  of  data  in  the  third  quadrant. 
However,  it  was  responsible  for  the  extension  of  the  5  C  contour  toward 
the  lower  left  corner  of  the  figure. 

The  isotach  analysis  in  Figure  6  shows  the  type  of  analysis  which 
can  be  generated  from  an  extremely  sparse  data  field.     The  most  signifi- 
cant point  in  this  analysis  is  the  manner  in  which  the  30  kt  contour  has 
been  distorted  in  response  to  the  36  kt  report. 
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The  stair-like  contour  across  the  bottom  of  Figure  7  was  caused  by 
the  lack  of  data  in  the  area.     The  field  was  initialized  to  O  C ,  and  since 
most  of  the  data  lay  in  the  vicinity  of  the  center,  none  of  the  points  lay 
close  enough  to  this  edge  to  be  spread  over  the  area  during  the  scanning 
portion  of  the  gridding  process. 

In  Figure  8  the  numerous  closed  centers  were  caused  by  the  separation 
of  data  points,  and  the  relatively  large  variation  in  values  between  these 
points.     This  850-mb  D-value  analysis  was  drawn  entirely  from  DROPSONDE 
data.     In  this  case  the  position  of  the  aircraft  at  the  time  that  the  drop 
was  made  has  a  great  influence  on  the  analysis,  and  is  probably  respon- 
sible for  the  fact  that  the  lowest  D-values  are  not  symmetric  about  the 
cyclone  center,  as  one  might  expect.     This  figure  also  exhibits  a  problem 
similar  to  that  shown  in  Figure  5  as  evidenced  by  the  closed  200-m  and 
100m  centers  located  approximately  60  nautical  miles  from  the  center  of 
the  cyclone.     In  this  case  the  problem  was  believed  to  have  been  caused 
by  the  relationship  of  the  position  of  the  densely  packed  negative  area  at 
the  center  of  the  figure  to  the  position  of  the  two  positive  65-m  reports 
(first  and  third  quadrants)  .     During  the  scanning  portion  of  the  gridding 
process  the  computer  generated  the  positive  centers  to  maintain  a  con- 
sistency in  the  gradient.     Here  again,  the  problem  appears  to  lay  in  the 
sparsity  of  the  data  and  should  be  aleviated  once  more  points  are  added 
to  the  data  set. 
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Figure  5.  3000-meter,  hurricane-stage  isodrosotherm  analysis  drawn 
from  23  data  points  using  a  5-C  contour  interval.  Scale:  1"=80  n  mi 
(1971  data). 
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Figure  6.     5600-meter,  tropical-storm  stage  isotach  analysis  drawn  from 
4  data  points  using  a  10-kt  contour  interval.     Scale:  lll=:80  n   mi     (1971 
data) . 
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Figure  7.  5600-meter,  tropical-storm  stage  isotherm  analysis  drawn 
from  11  data  points  using  a  5-C  contour  interval.  Scale:  1  "=80  n  mi 
(1971  data). 
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Figure  8.     850-mb,  hurricane-stage  D-value  analysis  drawn  from  18  data 
points  using  a  50-m  contour  interval.     Scale:  1"=80  n  mi       (1971  data). 
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VI.      CONCLUSIONS 

1.  The  computer  routine  developed  for  the  processing  and  analyzing  of 
aircraft  reconnaissance  data  performs  as  designed  and  is  now  ready  for 
more  extensive  use. 

2.  No  meaningful  determination  of  the  structure  of  eastern  tropical  North 
Pacific  (EASTROPAC)  cyclones  could  be  made  from  one  year's  data.     Never- 
theless,  1971  was  adequate  as  a  prototype  year  for  the  purpose  of  devel- 
oping the  computer  routine  for  processing  aircraft  reconnaissance  data. 

3.  The  majority  of  the  1971  data  collected  by  aircraft  on  cyclone  recon- 
naissance missions  in  EASTROPAC  were  taken  at  points  too  far  from  the 
center  of  the  vortex  to  be  of  use  in  a  study  concentrating  on  structural 
characteristics  of  the  relatively  small  EASTROPAC  cyclone,  as  evidenced 
by  the  small  number  of  points  found  to  be  valid  for  the  analyses. 
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VII.      SUGGESTED  TOPICS  FOR  FUTURE  RESEARCH 

1.  Since  the  computer  routine  appears  adequate  for  the  processing  and 
analysis  of  EASTROPAC  reconnaissance  data,  the  first  objective  of  future 
research  should  be  the  compilation  of  a  complete  data  set  using  all  recon- 
naissance data  gathered  from  1964  to  date.     (Such  is  already  in  progress 
as  a  follow-up  Master  of  Science  Thesis  in  the  Department  of  Meteorology, 
Naval  Postgraduate  School,  Monterey,  California.) 

2.  Once  a  complete  data  set  is  compiled,  analyses  at  (or  near)  each  of 
the  standard  pressure  levels  should  be  made  for  each  of  the  major  stages 
of  cyclone  development  (depression,   storm,  and  hurricane  stages). 

3  .      Satellite  interpreted  tropical-cyclone  data  and  ship-reported  data 
should  be  researched  and  used  to  increase  the  data  base. 
4.      Further,   in  connection  with  (1)  above,  the  natural  coordinate  system 
used  in  PROGRAM  ONE  might  be  extended  into  time  coordinates  as  well, 
using  as  a  base  value  the  time  at  which  the  particular  cyclone  reached  its 
maximum  intensity.     Each  reconnaissance  report  could  then  be  placed  at 
a  time  t  hours  before  or  after  maximum  cyclone  intensity.     This  would 
allow  more  precise  determination  of  the  structure  at  various  stages  of 
development. 
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